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We show that the Kijowski distribution for time of arrivals in the entire real line is the limiting 
distribution of the time of arrival distribution in a confining box as its length increases to infinity. 
The dynamics of the confined time of arrival eigenfunctions is also numerically investigated and 
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physical content of the Kijowsky distribution. 
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I. INTRODUCTION 

The problem of accomodating time as a quantum dy- 
namical observable has a long history and remains con- 
troversial to this day [1 S |1 H S i, II II H E3- The 
time of arrival (TOA) of a quantum structureless parti- 
cle at a point or a surface is one aspect of this problem- 
atic treatment of time in quantum mechanics in need of 
clarification 

Hll. Recently significant progress has been 
made, both at the operational and foundational fronts for 
the quantum time of arrival problem. Most notable at 
the operational level is the convergence of several analy- 
ses of the problem to the axiomatic distribution of arrival 
for freely moving particles due to Kijowski 6j. Kijowski's 
distribution was first obtained operationally (under some 
limiting conditions) by Allcock [3 and later shown to 
be arising from the quantization of the classical time-of- 
arrival @ ^| , inspirin g ex tension of the theory to the 
interacting case |0, IT^.Il4| or multiparticle systems . 
Most important of these developments in understanding 
the physical content of Kijowski distribution is the real- 
ization that, while axiomatic in nat ure, it can be obtained 
from an operational procedure [lfj, [lTJ , whose essence is 
to modify (filter) the initial state to counterbalance the 
disturbance introduced by the apparatus, in particular 
at low energies. 

At the foundational level, it has been shown that the 
non-self- adjointncss of the free time of arrival operator, 
widely regarded as due to the semiboundedness of the 
Hamiltonian in accordance with Pauli's theorem 0,13, 
can in fact be lifted by spatial confinement 0. Thus 
the concept of confined quantum time of arrival (CTOA) 
was introduced O . The CTOA-operators form a class of 
compact and self-adjoint operators canonically conjugate 
with their respective Hamiltonians in a closed subspace 
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of the system Hilbert space. Being compact, the CTOA- 
operators posses discrete spectrum, and a complete set 
of mutually orthogonal square integrablc eigenfunctions 
[2^|. The eigenfunctions of the CTOA-operators are 
found to be states that evolve to unitarily (i.e. according 
to Schrodinger's equation) arrive at the origin at their 
respective eigenvalues — that is the events of the centroid 
of the position distribution being at the origin and its 
width being minimum occur at the same instant of time. 

Now we are confronted with the problem of relating the 
already established results for the unconfined quantum 
particle to the confined one, in particular the question 
whether the established Kijowski time of arrival distri- 
bution is extractable or not from the confined time of 
arrival operators. Therefore, in this paper we give mean- 
ing to the limit of the discrete time of arrival distribution 
of the CTOA-operators defined on succesively larger seg- 
ments, proving that this limit is the Kijowski's distribu- 
tion. We begin by a short review of the properties of the 
confined and unconfined time of arrival operators; this is 
followed by numerical indications of the relation between 
these two cases, numerical indications that are strength- 
ened by the analytical proof of the following statement: 
The time of arrival operator on the full real line is the 
limit for I tending to infinity of the time of arrival op- 
erators defined for free motion in a segment of length 
21. The dynamics of the CTOA-eigenfunctions are then 
numerically investigated in the limit of arbitrarilly large 
confining lengths. We end by summarising and providing 
conclusions. 



II. THE QUANTUM FREE TIME OF ARRIVAL 
OPERATOR 

From a heuristic perspective, it is quite natural to as- 
sume that the quantum TOA may be associated with 
the quantization of the corresponding classical expres- 
sion. That is, if a classical free particle, of mass /i in one 



2 



dimension at initial location q with momentum p, will 
arrive, say, at the origin at the time T(q,p) — —[iqp -1 , 
then the quantum TOA-distribution must be derivable 
from a quantization of T(q,p), such as the symmetrized 



T = -^(qp 1 



P" X q) 



(1) 



Formally the time of arrival operator T is canonically 
conjugate to the free Hamiltonian, H = (2/i) _1 p 2 , i.e. 
[H,T] = ih. Equation- has been separately studied 
when the motion of the free particle takes place in the 
entire real line and when it is restricted to a segment of 
the real line. These studies led to some seemingly contra- 
dictory descriptions, in that the time of arrival operator 
is not self-adjoint for the full real line, while there is an 
infinite family of self-adjoint operators playing the corre- 
sponding role for the segment. Furthermore, for the full 
line case there is a well defined property of time covari- 
ance which is lacking in the confined case. 



A. The TOA-operator in H p 



L 2 (— oo, oo) 



In momentum representation, Eq. Q formally as- 
sumes the form 



T = 



ih/i 



2d_ 

p dp 



Under a sensible choice of domain [18|, T is a densely 
defined operator, and unbounded in H p :— L 2 (R, dp). 
It is however not self-adjoint, and admits no self-adjoint 
extension. It is in fact a maximally symmetric operator 
[l8j . According to the well-established theory for this 
kind of operator, there exists a decomposition of unity 
given by the following degenerate weak (or non-square 
integrable) eigenfunctions 



^(p) = Q(ap) 



ap 



Gvr//// I 



1/2 



a ip 2 1 j 2mh 



This is a family of functions parametrized by the eigen- 
value t and a discrete parameter a, which can be ei- 
ther + 1 or — 1, and gives the sign of the half- line on 
which the function presents its support in momentum 
space. It is straightforward to prove completeness, i.e., 

EJ^L^^tP'l^tp) = S(p-p'), and nonorthogo- 
nality, 



dp0(p)i>M(p) 



The physical content of the eigenfunctions is better ex- 
amined in coordinate representation [Tlj . Normalized 
wave packets with the form of quasi-eigenstates of T 
peaked at a given eigenvalue, but with a time width AT, 
have sharp space-time behaviour. Their average posi- 
tion travels with constant velocity to arrive at the origin 



at the nominal arrival time of the packet, in which the 
packet also attains its minimum spatial width. The pas- 
sage of probability density from one side of the origin to 
the other, is in summary as sharp as desired by taking 
AT — > As we shall see later, the unitary evolution 

of a (generalized) state leads to a distribution with 

point support at the origin at the instant corresponding 
to the eigenvalue t, thus lending support to its interpre- 
tation as a time of arrival eigenstate. 

Using completeness, we can now write the probability 
density for measured values of the T operator, 
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dp 
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which is, in fact, Kijowski's probability density [6|. 

An essential property of the distribution of arrivals in 
Kijowski's axiomatic approach is covariance under trans- 
formations generated by the Hamiltonian. It is evident 
that this property of covariance is indeed held by the dis- 
tribution above. Physically, it means that the probability 
of arriving at t for a given state is equal to the probabil- 
ity of arriving at t — r for the same state, once evolved a 
time t. This is the reflection on the probability density 
of the canonical commutation relation between H and T. 



B. The TOA-operator in Hi = L 2 [-l,l] 

In [ljj the naive definition in Eq. T for a free 
particle in a segment of length 21 is supplemented by 
adequate boundary conditions, leading to properly self- 
adjoint operators. Physically, the choice of time of ar- 
rival operator with spatial confinement of the particle 
in the interval [—1, 1} is dictated by the condition that 
the evolution of the system is generated by a purely ki- 
netic self-adjoint Hamiltonian, i.e. H = (2/i)~ 1 p 2 , where 
p is a self-adjoint momentum operator. This require- 
ment demands that the momentum operator p be one of 
the set {p 7 = —ihd q , \j\ < 7r/2}, with p 7 having the do- 
main consisting of absolutely continuous functions 4>{q) in 
Hi = [—1, 1] with square integrable first derivatives, which 
further satisfy the boundary condition <fi(—l) — e _2n '(f)(1) ■ 
Since T depends on the momentum operator, T is also 
required to be the corresponding element of the set of 
operators {T 7 }. 

In coordinate representation, T 7 becomes the Fred- 



holm integral operator (Tjip)(q) — f , T 7 (g, q') (p(q') dq' , 
for all tp(q) in H, where the kernel is given by 



T^ (q,q') 



"^fe^ (^H(q q') + e-*H{q' q)) , 



Ah sin 7 



(2) 
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T~,=o(q, q') 



(q + q')sgn(q - q) 



Aihl 



in which H is Heaviside's step function and sgn 
sign function. 

With this representation, one can show that 
and T 7 form a canonical pair 



(3) 
is the 



H 



of Tii — a non-dense 



over, the kernel T T (g. 

/V- 



subspace- 
q') of T„ 



7 

in a closed subspace 
every 7. 

i.e. 



-for every 7. More- 
1 7 is square integrablc 
\T 7 (q,q')\ 2 dqdq' < 00. This means that T 7 is 
compact, and, as a consequence, that it has a complete 
set of (square integrablc) cigenfunctions and its spectrum 
is discrete. This, it should be stressed, is a radically dif- 
ferent situation from that in the full line, where the op- 
erator has a continuous spectrum and is not self-adjoint. 

In what follows, we will need only the spectral prop- 
erties for the periodic confined quantum time of arrival 
operators (that is to say, 7 = 0). The operator T com- 
mutes with the parity operator. Furthermore, it changes 
sign under time inversion, which entails that its spec- 
trum is symmetric about 0. This suggests classifying 
its eigenfunctions in even and odd subspaces (which will 
be denoted by the subscripts e and o respectively), and, 
within each of those subspaces, by a discrete index n and 
the sign of the eigenvalue (indicated as a superscript). 
In this manner, and as computed elsewhere |l0| , the odd 
eigenfunctions are 



<o(<?) = Kqf 



s n q 
I 2 



with 



f ± (0 = e**e /4 [j_ 4 (0Ti/f(0 
while the even ones are given by 



with 



: (0 



e T^3/4 



J_#(0t*Ji(0 



Jdr n ) 



(4) 



(5) 



(6) 



(7) 



A n and B n are the normalization constants, while s n 
and r n are the solutions of the secular equations for the 
operator, namely the positive roots of J_ 3 (r) + |</| (r) + 

i Ji (r) = for the even case and of J_i/ A (s) = for the 
odd case. The eigenvalues are determined by 



rt = ± 



Ap n h ' 



(8) 



where p n stands for either r n or s„. 

It has been demonstrated in yjj that a CTOA eigen- 
function is a state that evolves unitarily — that is accord- 
ing to Schrodinger's equation — to being concentrated 



around the origin at its eigenvalue along a classical tra- 
jectory, i.e. a state in which the events of the centroid 
of the position distribution being at the origin and its 
width being minimum occur at the same instant of time 
equal to the eigenvalue. 



III. RECOVERING KIJOWSKI'S 
DISTRIBUTION FROM THE 
CTOA-OPERATORS 

A. Numerical examples 

We have thus seen that the time of arrival operators 
for confined and unconfined free motion have radically 
different properties. And it is not immediate how Ki- 
jowski's distribution follows from the spectral properties 
of the CTOA-operators. In this section we show that in- 
deed Kijowski's distribution is extractable from the dis- 
crete time of arrival distribution defined by the CTOA- 
operator. But before doing so we first provide further 
justification from numerical results. 

One basic difference between the two cases being ex- 
amined is that the spectrum is continuous for the full 
line, whereas it is discrete for the segment. As a con- 
sequence, the same can be predicated of the probability 
distributions for times of arrivals. Therefore, in order to 
compare like with like, it is convenient to use the accu- 
mulated probability of having arrived at the origin before 
a given instant for both cases. 

Also in order to perform sensible comparisons we shall 
consider initial states with compact support on the real 
line, which can also be defined for all segments which 
would include that compact support. In this manner the 
probability of having arrived at the origin before instant 
t can be computed for the same function both for the 
segment and the line. 

In the case of the segment, the probability of having 
arrived at the origin before instant t for the initial state 
■00 is calculated according to the standard quantum me- 
chanical prescription, 



(9) 



where y 7iS and t 7jS are the eigenfunctions and eigenval- 
ues of T 7 . This distribution is in turn investigated for 
increasing I, and compared with the accumulated proba- 
bility for Kijowski's distribution, 



IIk [t, ip ] dr. 



(10) 



where Uk [t, ipo] is the Kijowsky time of arrival density. 

As the length of the segment I increases, the discrete 
eigenvalues of the corresponding time of arrival opera- 
tors become denser, and we take advantage of this prop- 
erty to have another representation for comparison with 
Kijowski's distribution; namely plotting Kijowski's time 
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FIG. 1: The top-figure shows the accumulated probability of 
arrival versus time at the origin for a gaussian wave packet 
with mean momentum < P >= 100, full width at half max- 
imum a x = 0.05 and initial expected value < X >= —1 (all 
the quantities in atomic units). The dots correspond to con- 
fined motion with I = 10, and the solid line to the full line. 
In the lower figure we depict the corresponding probability 
densities (with the same notation), defined as explained in 
the text. 



of arrival density together with a discrete derivative of 
the probability of having arrived at the origin before a 
given instant, in the case of the segment, discrete deriva- 
tive that can be understood as a time of arrival density. 
This allows us an easier recognition of some features that 
would otherwise be obscured in the accumulated proba- 
bility. 

Figs, n 121 an d 01 show the accumulated probability 
for the discrete case and for the Kijowski distribution, 
together with their corresponding time of arrival densi- 
ties. In Fig. □ we depict a simple situation for an initial 
gaussian state, with good match between the distribu- 
tions for both the full line and the segment. Fig. [21 
corresponds to an initial state that presents the backflow 
effect |l9j . that is, that at some instants the quantum 
probability flux can become negative even if all the com- 
ponents of the state are of positive momentum. In the sit- 
uation depicted in Fig. [21 therefore, we can discriminate 
whether the discrete distribution associated with the seg- 
ment approaches Kijowski's distribution or, rather, the 
flux, which is also related to the density of arrivals. The 
result, as can be ascertained from the inset in Fig. [21 



FIG. 2: Accumulated probability of arrival (upper figure), 
versus time at the origin for two gaussian wave packet with 
mean momentum < Pi >= 200 and < P2 >= 100 with 
maximum interference at the origin, e.g., < Xi >= —1 and 
< X2 >= —0.5. The full width at half maximum is a x =0.05 
for both gaussians. The dots correspond to confined motion 
with / = 10, while the solid line corresponds to the full line. 
For completeness the dotted line depicts the quantum proba- 
bility flux at the origin, integrated up to the relevant instant. 
The inset shows the zone of maximal discrepancy between in- 
tegrated flux and accumulated Kijowski's probability, where 
the accumulated probability for confined motion matches Ki- 
jowski's. In the lower figure we show the corresponding proba- 
bility densities (with the same notation) defined as explained 
in the text, as well as the quantum probability flux at the 
origin: in this situation there is a backflow effect )l9|. 



is that Kijowksi's distribution is the one selected in the 
limit of large I. Fig. |31is a depiction of accumulated prob- 
abilities for different values of I, clearly showing conver- 
gence to Kijowski's accumulated probabilities; similarly 
with the discrete derivatives. 

All together, these numerical results suggest strongly 
that the I — » 00 limit of the time of arrival operators 
for the segment of length 21 tend to the time of arrival 
operator on the real line, in some sense. 



B. The limit for large I 

A first hint into the relationship between the confined 
and unconfmed TAO-operators is provided by the ker- 
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FIG. 3: The upper figure shows the accumulated probabil- 
ity of arrival at the origin versus time for two gaussian wave 
packets with mean momentum < Pi >= 5 and < P2 >= 1.5 
with < Xi >= -1 and < X 2 >= -0.5. The full width at half 
maximum is a x — 0.05 for both gaussians. The distinct lines 
correspond to different values of I: I = 3, long dashed line; 
I = 5, dotted-dashed line, and I = 15, circles. The solid line 
corresponds to Kijowski's distribution. In the lower figure we 
show the corresponding probability densities (with the same 
notation), as defined in the text. 



nel of the the later in position representation. In q- 
representation, the T operator takes on the form of the 
integral operator 



0»(<?) 



where the kernel is given by 



IT 



<p(q')dq' 



(q\T\q') = £ h (q + q')sgn(q 



(11) 



(12) 



Clearly this is the / — > 00 formal limit of the kernel 
To(<7, q'). It is also the limit as 7 — > of kernels T~(q, q'), 
and, since the large I limit washes out the effect of the 
boundary conditions, (i.e. the wavefunctions must nec- 
essarilly vanish at infinity) this implies that the kernels 
match up in the large I limit. For this reason, it will 
be sufficient for us to consider the periodic confined time 
of arrivals in the limit of arbitrarilly large confinement 
lengths. 



1. Kijowski distribution in position representation 

In what follows we will need the position represen- 
tation of the degenerate eigenfunctions of T. In q- 
representation, the eigenfunctions assume the form 



^f{q) 



,± iqP /H ( _p_y dp ^ (13) 

\2fj,n J 



where the integration is understood in the distributional 
sense. Even and odd combinations of these are likewise 
eigenfunctions with the same eigenvalue t, 



^ /o (g) = ^(^ + (g)±^-(ff)) 



(14) 



The full explicit expression for ip e J°(q) is obtained by 
direct substitution, leading to 



rt'°{q) 



,(1±2W 1 



iVVtXW exp 



■ m 

1 Aht 



(15) 



where D v is the parabolic cylinder function, and the even 
and odd case correspond to positive and negative sign 
respectively. 

In momentum representation it is known that these 
eigenstates form a resolution of the identity. It thus 
follows, owing from the unitarity of the Fourier trans- 
form, that the same happens in the position representa- 
tion Hilbert space, L 2 (R), i.e. f^^dttpf (q')ift (q) — 
6(q — q'), with j3 standing for e and o. They are likewise 
non-orthogonal. 

Because of the invariance of Kijowski's distribution un- 
der parity, it can be split in even and odd components of 
the wavefunction. In terms of tp\ \q)i it then assumes 
the form 



n^ (t) = n; o (i) + n^(i) 



where 



i Pt /o (q)i'o(q)dq 



(16) 



(17) 



In the following section, we will show that H^°(t), and, 
hence, Kijowski's distribution itself, can be obtained from 
the accumulated probability of arrival F$(t) in the limit 
I — > 00. 

We will do so by a proper identification of the limit 
of the eigenfunctions of the confined quantum time of 
arrivals for I approaching infinity as the eigenfunctions 
given by Eqs. (|15|l . To this end we will need the position 
representation of the eigenvalue problem for T, which is 
explicitly given in momentum representation in the form 



ijjih 1 



1 diPtjp) 
p dp 



tMp)- (18) 
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Multiplying both sides of the equation by p 2 and then 
Fourier transforming the resulting expression leads to the 
position representation of the above eigenvalue equation 



d 2 ip t (q) fiiqd(p t (q) 3/ii 



(19) 



where <fit(q) 



J: 



r ■ > - oo ex P ( 1 ? p) ^* W Straight- 
forward substitution of ip^°(q) in the differential equa- 
tion shows that they are linearly independent solutions 
ofEq. PI). 



and the operators T and To are invariant under parity. It 
follows that, if we rewrite an initial function in terms of 
an even and an odd components, their separate distribu- 
tions for times of arrival sum to the distribution for the 
total function, with no interference term being required. 
It is thus useful to separate the analysis in the even and 
odd sectors. That is, we need to examine whether the 
distributions F^ e ^°(t) (which are the separate contribu- 
tions to the probability of having arrived at the origin 
prior to instant t of the even and odd components of 
the initial state ?/>o) do indeed tend to f^^dt' U^°(t') 
respectively as I approaches infinity. 



2- F^(t) for large I 

rited out above, both 
full real line and in that of the segment the Hamiltonian given by 



First let us consider the contributions from the odd 
As we have pointed out above, both in the case of the eigenfunctions for the accumulated probability. This is 
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± ij-i 
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ipo{q) dq 
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where the sum includes only contributions from eigen- 
functions with eigenvalues smaller than t, I is large 
enough such that the interval [—1,1] contains the sup- 
port of tpo(q), and we have used the explicit value of the 
normalization factor A n . 

As our numerical computations demonstrate, for in- 
creasing / the dominant eigenfunctions contributing in 
the calculation of the accumulated probability come from 
large values of n. But as n increases the eigenvalues of the 
dominant contributors become denser. Therefore, for ev- 
ery time t and large I there exists a corresponding n(l,t) 
such that \t\ — T n (i t t) (that is, the difference between the 
eigenvalue closest to \t\ - for that value of I - and \t\) 
tends to as I tends to infinity. For any given r then we 
can write the corresponding s n as as \il 2 /4tH, with the 



adequate sign, in the inner integrand. 

As for the normalization factor, it is clearly the case 
that s n tends to infinity, in which limit we can use the 
asymptotic properties of Bessel functions and their zeroes 
to write l 3 J3(s n ) 2 sl/ 2 in the form 2l 3 /nsl/ 2 or, alterna- 
tively, (4/ 2 /7r) v /ft|r|/^. Finally, since the spacing be- 
tween the roots of J_i/4 tends to tt as n tends to infinity, 
we have the result that the spacing between successive 
values of r„ (with either sign) is AttHt 2 / (fil 2 ) , whence it 
follows that the sum X) T ± <t can ^ e substituted by the 
integral \xl 2 j l _ dr/AhnT 2 

Putting together all these asymptotic expressions, we 
see that Eq. (|2*U|l becomes 



p(i)o (t) _^ 2ft 



4 / 4 \4:Th 



tpa(q)dq 



(21) 



Differentiating this with respect to time yields the corresponding probability density 
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from which we extract the limit of the odd eigenfunctions 



as I approaches infinity. The limit is 

y?\q) 
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(23) 
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The relationship between eauations- ll5l and- l23l is not im- 
mediate. But their relationship can be established by 
substituting ip° dd (q) back in equation- fl^ , and finding 
that it is a solution to the differential equation, meaning 
it is an eigenfunction of the time of arrival operator in 
the entire real line with the eigenvalue t. Since ip° dd (q) is 
odd, then it must differ at most with y>^~ (q) by a constant 
factor. Expanding ip° dd (q) and tp^~(q) about q = 0, we 
find that they differ only only up to the irrelevant phase 
factor e~ l7r / 4 . Then we must have U° dd (t) = III (t) in 
the limit of infinite I. 

The same procedure can be carried out for the con- 
tribution of the even eigenfunctions, albeit with more 
cumbersome algebra, and the corresponding conclusion 
is derived. As a consequence, we have proved that the 
I — > oo limit of the discrete probability distribution for 
times of arrival for a spatially confined particle is indeed 
Kijowski's probability distribution. 

It might be thought that our proof is lacking in that 
above we assumed that the initial state has compact sup- 
port, while Kijowski's distribution must also be defined 
for other normalizable functions. However, our analy- 
sis can be extended to initial states with tails extending 
to infinity. Let ip(q) be such a state. There always ex- 
ists a sequence of states with compact supports, ip n (q), 
n = 1, 2, . . . , such that i> n (q) — * ip(o)- We can, say, pick 
ipi(q) and apply our above analysis to this initial state. 
Once we get the limit of infinite I, we follow it with the 
limit in n. The resulting limit is going to be Kijowski's 
distribution because of the continuity of the inner prod- 
uct. 



IV. THE DYNAMICS OF THE 
CTOA-EIGENFUNCTIONS IN THE LIMIT OF 
INFINITE I 

We have shown above that Kijowski distribution is the 
limit of the discrete distribution for arbitrarily large I. 
But what physical insight can we get from this realiza- 
tion? Recall that one of the surrounding issues against 
an ideal quantum time of arrival distribution is the fact 
that a quantum particle losses the localized property of 
its corresponding classical particle entity. Classically the 
concept of time of arrival is well-defined because a clas- 
sical particle has a well defined trajectory. This is con- 
trary to the fact that no such trajectory can be ascribed 
to the quantum particle — it has no definite position and 
momentum. 

Now the theory of confined quantum time of arrivals 
demonstrates that the quantum time of arrival problem, 
at the ideal level, i.e. at the level where measuring instru- 
ments play no explicit role, can be rephrased to finding 
states that unitarily arrive at a given point at a definite 
time — states in which the events of the centroid being 
at the origin and the position distribution width being 
minimum occur at the same instant of time. The QTOA- 
problem phrased in this way is well-defined because quan- 



0.04 

-TO- 

JZ 0.02 





1 2 3 4 5 

/(a.u.) 



FIG. 4: Width at half maximum (WHM) of the evolved odd 
and even eigenfunctions (upper and lower lines) at the corre- 
sponding eigenvalues closest to t = 0.01 versus length I. 



turn states have well-defined trajectories according to the 
Schrodinger equation. All these give us insight to the 
ideal physical content of Kijowski distribution. Let us 
see how. 

Let us consider some fixed time r, and for any given 
length I of spatial confinement, we can find an n such that 
the eigenvalue r„ is closest to r. Now consider a sequence 
of monotonically increasing Ps, 1%, 1%, I3, with l\ < I2 < 
I3 < .... Then there will be an n\ corresponding to l\ 
such that T ni is closest to r; and an n 2 corresponding to 
I2 such that t„ 2 is closest to r; and so on. For arbitrarily 
large lengths, we should have r ni ~ r„ 2 « .... « r, so 
that in the limit of infinite lengths they converge to r. 

We know that the eigenfunctions tp ni , ip n2 , ... will ob- 
tain their minimum variances at their respective eigen- 
values r ni , t„ 2 , .... How do the variances compare for the 
different </?„'s? Equivalently, what is the behavior of the 
variance with respect to tp n at the eigenvalue r« as I 
approaches infinity? In this case the variance increases 
with I. This is not surprising because the ip n 's would be 
thrown out of the Hilbert space for infinite I, i.e. they 
acquire infinite variances. 

However, if we substitute [width at half maximum 
(WHM) of the probability density \ip n (q,T n )\ at the 
eigenvalue r n ] for the variance, we find that WHM de- 
creases with increasing I. Figure-Q] demonstrate this. 
What is even more important is that the density 
IVn(<Zi T n)| a t the eigenvalue r„ tends, as I increases in- 
definitely, to a function with point support at the origin. 
Figure-El demonstrates this. These results imply that the 
CTOA-eigcnfunctions in the limit of infinite I (which are 
already outside of the Hilbert space) evolve to "collapse" 
at the origin at their respective eigenvalues. This allows 
us to interpret Kijowski's distribution as the ideal time 
of arrival distribution of collapsing states at the origin. 
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FIG. 5: Probability density |<£>n(i~)| 2 versus position at the 
corresponding closest eigenvalue r n to t — 0.01, for the even 
(upper figure) and odd (lower figure) eigenfunctions. The 
different lines are associated with I — 1 (thick solid line), 
I = 2 (dashed line), I = 3 (long-dashed line), I — 4 (dotted- 
dashedline) and I — 5 (thin solid line). 



V. DISCUSSION 

From a purely technical point of view, we have shown 
that Kijowski's (continuum) time-of-arrival distribution 
for free, unconfined quantum particles can be obtained as 
the limit of the (discrete) distribution that results from 
quantizing the time of arrival in a finite box, taking into 
account the eigenvalue spacing in the transition from 
sums to integrals. Note that the TOA operator is not 
self-adjoint in the unconfined case, but it is self-adjoint 
in the box. There are several possible readings of this 
fact. We will point out some of them, with the aim of 
opening a public discussion on the fundamental issues in- 
volved, rather than to settle final answers and exhaustive 
explanations. 

For those following a traditional, von Neumann's for- 
mulation of quantum mechanics, the connection with a 
self-adjoint operator may be satisfactory since, within 
this framework, all observables must be associated with 
self-adjoint operators. One of the customary roles of 
self-adjointness is to assure the orthogonality of eigen- 
functions so that, according to the projection postulate, 
repeated measurements would give the same result. It 
turns out however, that this idealization is not applicable 
to many observables and/or measurement operations. In 
particular, improper eigenfunctions for observables with 
continuum spectrum (such as momentum, energy or po- 
sition on the line) are not normalizable, and thus the 
collapse and repeated measurement idea cannot be ap- 
plied literally, even in principle, but only approximately. 
It seems then that, as long as some ideal probability dis- 
tribution may be computed unambiguously, there is no 



fundamental need to have a self-adjoint operator in these 
cases. Indeed, if one adopts the point of view that ob- 
servables are best formulated in terms of POVMs 0,0 , 
self-adjointness is not essential to describe observables. 
The eigenfunction orthogonality may however be seen 
anyway as a desirable, simplifying property, whereby a 
part (component) of the initial state is responsible for a 
given outcome (eigenvalue) and not for any other. What 
the present paper shows is that the time of arrival eigen- 
functions in the unbounded case differ, appart from a 
normalization factor, only outside the large box or far 
from the arrival point at the origin. The nonorthogo- 
nality in the continuum is thus caused by the distant 
eigenfunction behaviour with negligible overlap with the 
the initial wavepacket, and is thus physically irrelevant 
for computing the distribution. 

Another problematic matter is the interpretation of 
the discrete nature of the time of arrival eigenvalues in 
the box. A discretization of time variables should not 
be surprising. For any system with discrete energies or 
eigenmodes the corresponding time periods are also dis- 
crete. In the same vein, the TOA operator considered in- 
volves the discretized momentum operator in the denom- 
inator and therefore discrete eigenvalues. The problem 
thus is not in accepting the possibility of a discretization 
but in determining its operational meaning. "Measure- 
ments" , and "observables" are frequently highly idealized 
in quantum mechanics, to the point that all explicit refer- 
ence to an apparatus could disappear. This is useful, but 
may also let us without important physical references for 
its operational interpretation. Whereas an operational 
interpretation exists for Kijowski distribution in the con- 
tinuum, a direct operational interpretation of the discrete 
times of the confined case is still missing and is one of the 
challenges for future research. In any case, the discrete- 
continuum smooth transition found here may be a useful 
tool to generate new theories of, say, first time-of-arrival 
with interacting potentials, that can be later translated 
to the continuum and operationally interpreted or com- 
pared with existing operational proposals poT | . 
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